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We study the corrections to the galaxy three-point correlation function (3PCF) induced by weak 
lensing magnification due to the matter distribution along the line of sight. We consistently derive 
all the correction terms arising up to second order in perturbation theory and provide analytic 
expressions as well as order of magnitude estimates for their relative importance. The magnification 
contributions depend on the geometry of the projected triangle on the sky plane, and scale with 
different powers of the number count slope and redshift of the galaxy sample considered. We evaluate 
all terms numerically and show that, depending on the triangle configuration as well as the galaxy 
sample considered, weak lensing can in general significantly contribute to and alter the three-point 
correlation function observed through galaxy and quasar catalogs. 
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I. INTRODUCTION 

The detection and measurement of the three-point and 
higher order correlation functions in position space and of 
their Fourier counterparts in the large-scale galaxy distri- 
bution date back to the late seventies and early eighthies, 
see, e.g. [1-4]. At that time such observations served 
in the first place as a confirmation of the hierarchical 
model of structure formation. Cosmological perturba- 
tion theory, for instance, makes specific predictions for 
matter higher order correlators which are generated by 
non-linearities in the equations describing gravitational 
evolution, even when Gaussian initial conditions are as- 
sumed [5, 6]. An extra source of non-Gaussianity is given 
by non-linearities in the galaxy bias, i.e., in the relation 
between the observed galaxy distribution and the under- 
lying dark matter fluctuations [7]. 

The measurement of the galaxy bispcctrum or of the 
three-point function in redshift surveys can provide a de- 
termination of linear galaxy bias independently of the 
normalization of the dark matter perturbations [8~10], 
and, in the framework of the halo model, on the parame- 
ters controlling the halo occupation distribution [11, 12]. 
This kind of analysis has been performed for current 
large-scale galaxy surveys such as the 2dF Galaxy Red- 
shift Survey (2dFGRS) [13-15], and the Sloan Digital 
Sky Survey (SDSS), [16, 17]. More recently, it has been 
shown that the bispcctrum can also provide insights into 
the non- linear nature of galaxy bias which can be used to 
constrain models of the power spectrum in the range of 
the acoustic oscillations [18-20]. In general, it has been 
shown that in surveys such as the SDSS, adding bispec- 
trum information allows for an improvement in the de- 
termination of cosmological parameters as compared to 
using the power spectrum alone [21]. 

An additional contribution to the galaxy three-point 



function is expected from primordial non-Gaussianity, i.e. 
if perturbations at early times arc not perfectly Gaus- 
sian [22-28]. Current constraints on primordial non- 
Gaussianities derived from CMB observations, [29-33] 
still allow for a component potentially observable in fu- 
ture large-scale galaxy surveys, [27, 34]. 

Recent interest in the observation of acoustic features 
in the galaxy correlation as a means to accurately map 
the expansion history of the Universe and thereby con- 
strain dark energy parameters [35-38], has motivated 
several proposals for future large surveys at high red- 
shift, as, for instance, the Hobby-Eberly Dark Energy Ex- 
periment (HETDEX) [39], the Wide-Field Multi-Object 
Spectrograph (WFMOS) [40], the Advanced Dark En- 
ergy Physics Telescope (ADEPT) mission, the ESA Eu- 
clid mission [41], or the Cosmic Inflation Probe (CIP) 
[42]. 

One of the technical and theoretical challenges that 
these projects will have to face is that the galaxy power 
spectrum or correlation function of high-redshift objects 
is affected by weak lensing magnification due to the in- 
tervening matter along the line-of-sight [43-47]. When 
a magnitude limited survey is carried out, weak gravita- 
tional lensing acts on what is observed in several ways. 
The first effect arises because a congruence of photons 
traveling from a distant galaxy to the observer is focused 
or defocused by gravitational potentials along the way, 
thus resulting in magnification or demagnification of the 
object observed. Weak lensing therefore pushes objects 
over the threshold magnitude of the survey, in cither di- 
rection. The resulting catalog does not contain all the 
objects that are intrinsically brighter than the limiting 
magnitude, and contains some objects that arc intrinsi- 
cally fainter than this magnitude. The size of this effect 
depends on the slope of the magnitude distribution of 
galaxies: the steeper the slope, the stronger the effect of 
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weak Icnsing magnification on the observed number den- 
sity will be. The second effect is due to the fact that weak 
lensing also stretches or contracts patches of the sky, thus 
varying the density of objects observed. Through these 
two first-order effects, lensing affects the observation of 
galaxy overdensities. Since these are used as a way of 
quantifying matter overdensities, weak lensing enters the 
measurement of any correlation function that is based on 
a magnitude limited survey. We do not consider the ef- 
fects of the lensing-induced displacement of sources here 
(i.e. the "smoothing effect" ), as we expect it to be smaller 
than the magnification contribution. 

This work attempts a first step towards quantifying 
and possibly exploiting the weak lensing contribution to 
the three point correlation function of pointlikc objects 
like galaxies or quasars. 

The intrinsic galaxy three-point function depends on 
three parameters (e.g., the three sides of a triangle). The 
lensing effects induce an orientation dependence in the 
observed 3PCF, specified e.g. by two angles. Hence, the 
lensing effects depend on five parameters, in addition to 
the number count slope and the rcdshift considered. In 
this paper, we do not attempt to fully explore this vast 
parameter space. Rather, after presenting the full set of 
equations, we focus on specific parameter sets chosen to 
allow for an understanding of the basic features of the 
lensing effects. 

For this work, we use the framework of second order 
cosmological perturbation theory [5, 6]. While the valid- 
ity of this approach is limited to large scales, it allows 
for a straightforward and consistent analytical treatment 
of the magnification correction. Moreover, the results 
are physically intuitive, and, as we will show, all relevant 
expressions reduce to manageable integrals over a small 
number of auxiliary functions. 

The paper is organized as follows. In Sec. II we present 
the un-lensed 3PCF, we then derive all the terms arising 
because of weak lensing and provide order of magnitude 
estimates for the different correction terms. In Sec. Ill we 
present results for four specific triangular configurations 
and for the reduced three-point function. We also obtain 
an estimate of the observability of the lensing effect for 
future surveys. We conclude in Sec. IV. The details of 
the calculation as well as discussions of approximations 
used are relegated to the Appendices. 



II. MAGNIFICATION OF THE 3PCF 

A. The un-lensed galaxy 3-point correlation 
function 

We want to study the effects of lensing magnifica- 
tion on the observable galaxy 3-point correlation function 
(3PCF) on large scales. A non- vanishing three-point cor- 
relation function, even for Gaussian initial conditions, is 
expected due to non-linearities in the dark matter evolu- 
tion, as well as due to non-linearities in the galaxy bias. 



i.e. the relation between the observed galaxy distribution 
and the underlying dark matter field. 

Here we consider the expression for the dark matter 
3PCF predicted by tree-level perturbation theory (PT) 
[5], for a review see [6]. The dark matter bispectrum, 
with Gaussian initial conditions, is then given by 

B(fci,A:2,fc3) = F2(ki,k2)PL(fci)PL(fc2) + 2pcrm., (1) 

where PL{k) is the dark matter linear power spectrum 
and the PT kernel F2 is given by 

The three-point matter correlation function in position 
space is defined as 

C(X1,X2,X3) EE (,5(X1)(5(X2)(5(X3)) 

= Jd^kie-'^' '''j'd^k2e-'^' ''''B{ki,k2). (3) 
Here, xi, X2, X3 denote the vertices of the triangle; 

^ij — : '^ij — \ 1 1 (^) 

and (5D(ki23) is the Dirac delta fimction for the sum of 
the three wavevectors ki23 = ki + k2 -I- k3 . 

A simple expression for X2, X3) has been obtained 
by Jing and Borner [48] assuming the tree-level matter 
bispectrum given in equation (1), 

C = ye(^i2)C(xi3) 

- ['72 (2^12 )f?o (3^13) + ?7o(a;i2)'72(a;i3)]xi2 • X13 

+ y [e(xi2)e(xi3)(xi2 •xi3)^ 

+'y2(a;i3)e(xi2)xi2 + ?y2(a;i2)e(a;i3)a;i3 

+3772 (2:12 )?72 (3:13)] + 2 perm., (5) 

where the auxiliary functions 77; (x) and e{x) are defined 
in Appendix A. They are related to derivatives of the 
two-point correlation function. 

The matter 3PCF is usually measured using survey 
data, under the assumption that galaxy overdensities 
trace matter overdensities. A simple, local model for the 
galaxy bias, consistent with the perturbative treatment 
of dark matter evolution, is given by the expansion in 
terms of the matter overdensity J(x) [7], relating (up to 
second order) the galaxy overdensity (5g(x) to 5(x) as 

J,(x)^6i<5(x) + |52(^)^ (g) 

where the two constants bi and 62 are the linear and 
quadratic bias parameters, respectively. In order to be 
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consistent with the PT expression for the dark matter 
bispectrum it is necessary to retain contributions propor- 
tional to 62. This is because the galaxy 3PCF receives 
two contributions of the same order: the first one from 
the dark matter, tree-level, 3PCF and the second one 
from the non- linear bias parameter 62- From equation (6) 
one obtains 



where k(x) and 7(x) represent, respectively, the conver- 
gence and shear (see equations (32) and (33) and ap- 
pendix B for the definitions). Following the perturbative 
approach adopted for the estimation of the large-scale 
3-point correlation function, and using the fact that k 
and 7 are linear in the matter density 6, we expand the 
magnification factor in equation (11) to second order as 



bl C(X1,X2,X3) 

blb2 [62 a3 + 2pcrm.], (7) 



fi 1 + 2k + 3k^ 



(12) 



where C(xi , X2, X3) is the matter 3PCF and ^ij is a short- 
hand notation for the matter two-point function for a 
given separation 



From equations (9, 12), the observed galaxy number den- 
sity can then be expressed as 



= ^(x„) ((5(xO(5(xj)). 



B. Magnification Corrections 



(8) 



Due to the effect of weak lensing magnification, the ob- 
served galaxy number density n{< m,n) below a certain 
limiting magnitude m in a given direction n in the sky 
is not equal to the actual density ng, but is related to it 
via [49] 



n(< m, n) ~ n(< m) ( 1 + fei (5 + 



X [I + Ci K + ^C2 + ^ci -f"^) , (13) 



where we have defined the constants ci = (5s — 2) and 
C2 = ci(ci + 1) = (5s — 2)(5s — 1). Up to second order in 
(5, the observed galaxy overdensity is then 



2.5S-1 



(9) 



where no = n[l -I- dg{n)], fi{h) is the magnification, and 
s is the slope of the number density as a function of the 
magnitude m 



5ofcs(x) ~ bi S + ci K + -62 5^ + bi d S K 



-^C2K 



2^17 



(14) 



_ dln[no(m)] 
dm 

The magnification ^ is given by 



1 



(1-^)2-72^ 



(10) 



(11) 



From this expression one derives sixteen different con- 
tributions to the observed 3PCF, all corresponding to 
terms of the same order in perturbation theory. Making 
use of the short-hand notation 5i = '5(xi), = K(xi) and 
7, EE 7(xi), we have 



((5ohs(xi)Jofcs(x2)(5o6s(x3)) = bl ((5i(52(53) + b2 [ - {5i525l) + 2 perm.] 

-]-bl ci [((5i(52K3) + 2 perm.] + bi 62 ci \^{h5lK3) + 5 perm.] 
+bl ci [((5i(S2(53K3) + 2 perm.] -f bl C2 \^{5i52nl) + 2 perm.] 
+^1 ci [^('5i'5273) + 2 perm.] + bi cj [{S1K2K3) + 2 perm.] 
+^2 cj [^((5?K2K3) -I- 2 perm.] + b\ c\ [{5i52H2Hz) + 5 perm.] 
+61 ci C2 \^{5iH2nl) + 5 perm.] + bi c\ \^{5iH2ll) + 5 perm.] 
+Ci (K1K2K3) + c\ C2 [^(«iK2K3) + 2 perm.] 

+Ci \^{f^ii^2ll) + 2 perm.] -I- 61 c\ [(^iKiK2K3) + 2 perm.]. (15) 
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Each of the contributions may involve several permu- 
tations of vertices. Note that, as one takes into ac- 
count second order terms in the convergence, one can- 
not distinguish unambiguously between "galaxy-galaxy- 
lens", "galaxy-lens-lens", and "lens-lens-lens" contribu- 



tions. For practical purposes only, we will nevertheless 
identify four groups of terms, the first, referred to as the 
"galaxy-galaxy-galaxy" group (or GGG), being given by 
the intrinsic matter 3PCF plus contributions due to non- 
linear bias: 



bl {S1S2S3) = bl C(xi,X2,X3) 



[GGG 



bl 62 2 ('^i'^2^I) + 2 perm. 



— bi &2 (C12 ^13 + C23 + Ci3 $23) [GGG — B]. 



(16) 
(17) 



For the rest of the paper, we assume without loss of 
generality that xi ^ X2 ^ XSi where Xi represents the 
comoving distance to the point x^. In the context of 
the Limber approximation, we can then neglect all terms 
where lenses in the direction of closer vertices are cor- 
related with overdensities further away, such as (ki(S2), 



(kiK2(53), and so on. The applicability and limitations of 
the Limber approximation in this context are discussed 
in Appendix C. 

We can now express the non-vanishing contributions 
due to lensing magnification as follows, distinguishing a 
"galaxy-galaxy-lens" (GGL) group. 



b\ ci [{5i52H3) ^ 


f 2 perm. 


] = bl 


Cl 


(<5i(52Ai3) [GGL -A], 






(18) 


2 Ci 


2('5l^2'^3)^ 


- 5 perm. 


= 61 


62 


Cl [^12(('52'«3) + {SiK^)) - 


^Cl3('5lK2)] 


[GGL-B], 


(19) 


bl ci 


[((5i(52^3K3) ' 


f 2 perm. 


1 = bl 


Cl 


(?23(^lK3) + 63(^2K3) + 


i2z{5lH2)) 


[GGL-C], 


(20) 


1 C2 




- 2 perm. 


= bl 


C2 


(<5lAJ3)(^2A^3) [GGL- 


D], 




(21) 


1 Cl 




- 2 perm. 


= bl 


Cl 


(<5i73)('5273) [GGL- 


E], 




(22) 



a "galaxy-lens-lens" (GLL) group. 



bl cl [((5iK2K3) + 2 perm.] = bi cl {S1K2K3) [GLL — A], 



b2 cl 



{S1K2K3) + 2 perm. 



bl q [{6162K2K3) + 5 perm. 



bl Cl C2 

bl cl 



-{6iK2nl) + 5 perm. 
^('5iK273) + 5 perm. 



b2cl {dlK2){dlK3) [GLL-B], 

bl cl [{5lK2){S2K3) + (C12 +Ci3)(k2K3) 

+ (62 + 63)(kiK3) + (^13 + 63)(kiK2)] [GLL - C], 

bl Cl C2 [(((5lK3) + {SlK2)){n2K3) + {S2H3) {K1K3)] [GLL - D], 



and a "lens-lens-lens" (LLL) group. 



C? C2 



cl 



q {K1K2K3) 

K1K2K3) 4- 2 perm. 
K1K27I) + 2 perm. 



bl cl [{S1K1K2K3) + 2 perm.] 



cl (K1K2K3) [LLL - A], 

Cl C2 {{kiK2){k.iK3) + («;iK2)(k2'«3) + (K1K3) (K2K3)) [LLL - B], 

cl i{'JlK2){llK3) + (ki72)(72K3) + (ki73)(k273)) [LLL - C], 

bl cl (((5iK2)(kiK3) + (<^iK3)(kiK2) + {S2K3){kiK2)) [LLL - D]. 



(23) 
(24) 

(25) 
(26) 



= bl cl [((5i73)(«:273) + (-^i 72) (72^3) + ('5273) ('«i73)] [GLL - E], (27) 



(28) 
(29) 

(30) 
(31) 
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The right hand sides in the equations above include 
exphcitly aU non-vanishing permutations. Weak lensing 
and non-Unear bias each couple terms which, as "dis- 
connected" terms, would not contribute to the observed 
3PCF otherwise [6]. The non-linear bias couples the mat- 
ter overdensity at a fixed location with itself, see e.g. the 
GGG-B contribution [equation (17)]. The magnification 
effect couples a lens along the line of sight to a galaxy 
with that galaxy, with a "coupling strength" given by 
ci = 5s — 2. For example, for the first term of the GGL- 
C contribution [equation (20)], a galaxy at xi correlating 
with a lens along the line of sight to the galaxy at X3 , and 
a galaxy at X2 correlating with the one at X3 , result in a 
contribution to the observed galaxy three-point function. 
In addition, there are mixed terms where one coupling is 
due to non-linear bias while the other is due to lensing 
[e.g., GGL-B, equation (19)]. 



C. Evaluation of Corrections 

We now proceed to evaluate the terms contributing 
to the observed 3PCF, equations (16-31). As shown in 
Appendix B, the magnification and shear perturbations 
are given as the following line-of-sight integrals: 



«;(x,) = ::a„i?o' / \-^^^^t¥^'5 (y.(x);x) , (32) 



7(xj) 



a(x) 



2 l""' ^^WLiX^,X) 



a{x) 

3 fej - + 2ikik2 



'5k(x)e 



-«k-yi(x) 



(33) 



Here, WLixi,x) = xlXi{x% - x), and Xi is the distance 
to Xi, while X is the lens distance. yi(x) = (x/xO^j 
denotes a point moving along the line of sight connecting 
the observer and the source located at x,; at a distance 
Xi- In the following, we also use the shorthand notation 

C = SflmliQ /2. 

In addition to the matter two-point and three-point 
functions ^ and C, we need to calculate permutations 



of 3 different connected three-point terms, (SSk), (Skk), 
(kkk)^ and of four different two-point terms, (Sk), (Sj), 
(kk), (kj)- It is quite remarkable that they can all be 
calculated with the help of only four auxiliary functions, 
as detailed in Appendix A. 

We will make use of the Limber approximation [50] in 
the evaluation of the magnification terms. Each appli- 
cation of the approximation sets the lens distance x in 
equations (32) and (33) to the value of the observable 
(e.g., 5) correlated with, thereby reducing the number of 
line-of-sight integrals by one, and fixing the separation 
r to r±, the projected (or transverse) separation. See 
Appendix C for a discussion. 

By defining the "transverse" two- and three-point cor- 
relation functions, ^ and C, respectively, and the auxil- 
iary function e (see Appendix A for more details), one 
obtains very similar expressions for almost all magnifica- 
tion terms. The projected functions are given by 



|(xi,X2;x) 
C(xi,X2,X3;x) 



xPL{kr,x), (34) 
{2tt)^J d^ki^e-'^'^-^'''^-'''^^ 

X J (i2fc2±e-*^-L-(''^^-^^-L) 
xB(kii,k2i;x), (35) 



and 



e{x;x) 



(27r) 



dkPLik;x)[kxJo{kx)~2Ji{kr)]. (36) 



We give all relevant quantities in terms of these functions 
below, first as the exact integral, then using the Limber 
approximation (L.A.). In all instances, Xij± denotes the 
transverse component of the separation, i.e. the projec- 
tion onto the sky plane. 



((5(xi)(5(x2)k(x3)) 



/^j WLiX3,X)./ , ^^ 

<^ / <^X —X C(xi,X2,y3(x)) 





L.A. ^Wl(x3,Xi 



C 



a(x) 

C(xi±,X2±,X3i;Xi) (5n(xi - X2) 



a(Xi) 

2 i""^ , Wl{x2,x) r , ,Wl{x3,x') 



((5(xi)a^(x2)k(x3)) ^ I dx , , 

Ja a(x) 

^2^l{X2,Xi)Wl{x3,Xi 



dx' 



a^(xi) 



a(x'y 

C(xix,X2±,X3_l;xi) 



C(xi,y2(x),y3(x')) 



(37) 



(38) 
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.3 r ^ WL{xi,x) WLix2,x )WUx3,x)' 
a^x) 



I dx ^^^^'^^ '^''^''^' C(yi±(x),y2±(x),y3±(x);x), (39) 



(<5(xi)«(x2)) = C^'^'dx ^^jg^^ ^(xi,y2(x)) 

=i C -, — ^ ^ Xix,X2±;xi), (40) 

a(Xi) 

/ . w ^\ ^ ^£(Xi,x) , W^l(X2,xO ^, . x , 

(k xi)k xa) = C / dx 1—. / o^X i—R ? yi X 'y2(x)) 

Jo a(x) Jo a(x') 

^2£^ ^x ^^^"^'^"i(^^^"""^ e(yi.(x),y2.(x);x), (4i) 

(^(xi)7(x2)) = ^^^^^^ Jd^kPUKx f^ - % ^-^^^^ e-My.(x)-x.) 

LA. Wl(x2,Xi) 2 N /.o^ 

- C* — ^ — a;i2,±e(xi2,±;xi), (42) 

a(Xi) 

/ , w ^\ r<2 T'^ W^l(xi,x) f''\ ,Wl{x2,x') 
(k(xi)7(x2)) C / dx 1—. / dx 



a(x) Jo a(x') 

rminixux.) WUXU X)Wl{X2, x) 2 -r . . ^ .."^^ 

/ c^x j^r^^ yi2._Le(yi2,x(x);x)- (43) 

Jo a (x) 



r 



A few comments are in order. To gain a better un- 
derstanding of the effect of applying the Limber approx- 
imation, consider equation (38). The first, exact, line 
requires the integration of the matter three-point func- 
tion - weighted by the appropriate geometrical lensing 
factors - over all triangles that have one vertex fixed at 
xi while the other two (located at y2 and ya) are mov- 
ing along the lines of sight connecting the observer to the 
source locations at X2 and X3, respectively. Similarly, in 
equation (39) we integrate the unlensed 3PCF over all 
triangles with vertices (yi,y2,y3) lying on the three lines 
of sight connecting the observer to the sources. For all 
these terms, the Limber approximation counts the largest 
contributions to the integrals, which arise from configura- 
tions that are perpendicular to the line of sight, i.e. where 
the separation Xj^n along the line of sight is of order x^ j^ 
or less. This is because the two and three point corre- 
lation functions increase as the separations get smaller. 
The largest contributions to the integrals will therefore 
arise when the distance between the points y, moving 
along the lines of sight are smallest, which corresponds 
to transverse separation vectors, i.e. perpendicular to 
the line of sight. Hence, after the Limber approximation 
all terms are proportional to the "transverse" correlation 
functions evaluated at the transverse displacements x^j^, 
yn_ . Clearly, the magnification terms depend on the size 
and configuration of the projected triangle, rather than 



the three-dimensional triangle itself. 



However, it is important to note that the Limber ap- 
proximation has limitations. Most severely, for the GGL- 
A term [equation (37)], applying the Limber approxima- 
tion yields an unphysical (5d(xi ~X2), as the lens is forced 
to be at the same distance x as both 5i and 52- To avoid 
this artifact of the approximation, we evaluate the exact 
integral for this term. For all other terms, wc discuss the 
accuracy of the Limber approximation in Appendix C. 
Briefly, we found that, whenever an integral over the line 
of sight remains in the final expression, the Limber ap- 
proximation is accurate. This is the case for the (kkk), 
(kk), and (nj) terms. For the other terms, the approxi- 
mation can become invalid in certain limits. 



After specifying the vertices of the triangle, equa- 
tions (37-43) above then allow a calculation of all the 
terms contributing to the observed galaxy 3PCF. The 
results of this calculation are reported in Sec. HL Be- 
fore plunging into the full calculation of all the terms in 
equations (16-31), however, it is useful to assess the rel- 
ative magnitude of the different terms, in order to have a 
better understanding of which of them provide the most 
relevant contribution to the observed 3PCF. 
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D. Order of magnitude estimates 

We will give an estimate of the order of magnitude of 
the various contributions in terms of the matter corre- 
lation function ^(r), or, equivalently, the dimensionless 
matter power spectrum A^(fc), since ^{r) ~ A^(l/r) (see 
also section 3 of [44]). 

In general, we consider a single scale r for the trian- 
gular configuration, so that we might expect these esti- 
mates to be more accurate for equilateral configurations. 
We let r± denote the typical extent of the triangle in the 
direction transverse to the line-of-sight, while r|| repre- 
sents the longitudinal extent. The relative magnitude of 
these is related to the orientation of the triangle, with 
rj^ = r, r|| =0 for a triangle in the sky plane, and 
rj_ ~ 0, r|| ~ r for a triangle perpendicular to the sky 
plane. 

We approximate the matter 3PCF simply by the hier- 
archical relation ({r,r,r) ~ £,'^{r) so that the (SSS) con- 
tribution is given by 

C,g, ^ bliS^S^Ss) ^ bieir). (44) 

The contributions due to magnification involve the pro- 
jected functions ^ and ( which can be approximated as 

ar±) - r^e(r±) r^A^l/r^), (45) 

and 

ar±)^rie{r±)^rl[A\l/r^)r. (46) 

Note that ^ and C, as opposed to ^ and C: are not di- 
mensionless, but rather have dimensions of length and 
length squared, respectively. We use these approxima- 
tions in equations (37-43). Whenever there is an integral, 
we pull the correlation functions out, setting x ~* Xi/2, 



yi± r±/2, and integrate over the lensing weight func- 
tions (which can give appreciable numerical factors). 
When not integrated over, the lensing weight function 
is approximated as 

WL{x2,Xi)~X2-Xi=r\\. (47) 

As an example of this approximation scheme, consider 
the GLL-A term in equation (23), 

OlCi(diK2K3) = OlCiC — C 

a^iXi) 

- hclC^il + zfrlrie{r±), (48) 

where the first equality follows from the Limber equation 
in equation (38) and the second line implements the order 
of magnitude estimates. Here, z denotes the redshift at 
which the triangle is located, and the correlation func- 
tions arc to be evaluated at that redshift. Recall that 
C oc Hq and the canonical GGG term is of order 6f (r) , 
so the ratio of this contribution to the canonical term is 

GLL-A /ci, .C^{r±)V , . 

^^-(^(l + ^)(i?on,)(i^or.)^j • (49) 

The order-of-magnitude expressions can be used to or- 
der the contributions according to significance. Since 
HqX is of order 1, the main suppressing factors here are 
terms like Hi^r^ while the ratios of correlation functions 
can enhance the expressions somewhat (since r ^ < r and 
^ rises rapidly for small r). The GLL-A contribution is 
suppressed by a factor of order (r/x)^, denoted 0{^) 
below, and so is expected to be very small. 

Applying this approximation to all contributions leads 

to 
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[LLL-C 0(2)], 

[LLL-D 0(5)]. 



The order in powers of Hgr is indicated by the symbols 
0(71) in each hne. Comparing with the exact expressions, 
we have verified that the above expressions correctly re- 
produce the relative ordering of the different terms. The 
largest contributions to the magnification corrections of 
the 3PCF arise from the 0(1) and 0(2) terms. Terms 
of order 0(5) and 0(4) arc significantly suppressed. 

The significant GGL and GLL terms do not increase 
dramatically with redshift, as they scale with (1 + z) or 
(1 + z)2. The terms belonging to the LLL group on the 
other hand scale with third and fourth power of the red- 
shift: one expects their weight relative to the other terms 
to increase with increasing redshift. This is consistent 
with a similar behavior of the GL and LL magnification 
terms in the two-point correlation function, as already 
noticed by [43, 44, 47]. 



III. RESULTS 

The relative magnitude of the magnification correc- 
tions to the galaxy 3PCF derived in the previous section 
is highly dependent on the shape of the specific triangu- 
lar configuration considered, as well as on its orientation 
with respect to the line-of-sight. Typically, as shown by 
the estimates presented above, we expect a large effect for 
triangles whose projection onto the sky is much smaller 
than the three-dimensional extent, i.e. elongated trian- 
gles oriented along the line-of-sight. Another factor to be 
taken into account is the specific value for the linear and 
quadratic bias parameters as well as the number count 
slope that we can expect for the sample under consider- 
ation. 



(57) 
(58) 
(59) 
(60) 
(61) 
(62) 
(63) 

I 

A. Specific triangle configurations 

We will not present here a complete analysis of the de- 
tectability of this effect as this is beyond the scope of the 
present work. We will limit ourselves to consider some 
specific configurations at different redshifts and compute 
the relative magnification contributions. In section III C 
we will comment on the uncertainty of the 3PCF that 
one might expect in future high-redshift surveys for these 
specific triangles. 

We consider four classes of triangular configurations. 
The first three correspond to triangles lying in a plane 
that contains the line-of-sight. In this case, the projec- 
tions of the three galaxy positions on the sky plane will 
be aligned. We fix the transverse separation between 
the points and vary in different ways their line-of-sight 
separation. Let us use r^^x to denote the transverse sep- 
aration between galaxies at positions and x^, while 
Tyji = Xj ~ Xi is the separation along the line of sight. 

We carefully chose triangle configurations and scales to 
ensure that the trec-lcvel perturbation theory approach 
remains valid. We can assume loop corrections in per- 
turbation theory to be subdominant above scales of the 
order of 20ft.^^Mpc. At the same time, while the lous- 
ing signal becomes relatively larger at large scales, it will 
still be difficult to measure the three-point correlation at 
scales of hundreds of Mpc, even in forthcoming redshift 
surveys. 

For the first triangle configuration ("Tl"), schemat- 
ically represented in the left panel of Fig. 1, we keep 
ri3,i_ = r-23,i = 20 h^^ Mpc, ri2,± = 40 /i^^ Mpc and 
X2 ~ Xi = 20/i~^Mpc and we evaluate the 3PCF and 
its magnification corrections as the line-of-sight separa- 
tion r = X3 ~ Xi of the third galaxy is varied from 20 to 
200/1-1 Mpc. 
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FIG. 1: Triangular configurations Tl and T2. The dashed 
(red) line shows a second triangle of the same class. 



Ill the second class of configurations ("T2") we keep 
the same values for the transverse separations as for Tl, 
but we consider the position X2 to be halfway between 
Xi and X3 (see Fig. 1, right panel), with Tiz.\\ again vary- 
ing between 20 and 200/1"^ Mpc. The third one ("T3") 
is analogous to Tl but with fixed r23j| = 20 /i^^ Mpc 
and ?'i2.||, ''is, II varying. All these classes represent a 
continuous set of triangles starting from an overall size 
of tens of mcgaparsec which become progressively elon- 
gated and aligned with the line-of-sight as r\\ increases. 
From equations (37-43) above, we expect an increasing 
contribution from magnification for such elongated con- 
figurations: while the intrinsic (unlcnsed) correlation de- 
creases with the triangle size, the Icnsing terms will re- 
main roughly constant, since they depend on the 'pro- 
jected triangle which we keep fixed. 

In addition to these three peculiar classes of configu- 
rations, we evaluate the lensing contributions for a set of 
equilateral triangles of increasing size, lying on a plane 
subtending a 45° angle with the line-of-sight (see Fig. 2). 
One can take this to be an example of a triangle with 
generic orientation, where the parallel and perpendicular 
separations increase proportionally. The relative lensing 
contributions are in this case expected to be nearly inde- 
pendent of the side r of the triangle. 

Assuming a flat ACDM cosmology with parameters: 
a™ = 0.27, VLb = 0.046, h = 0.72, n, = 1, ag = 0.9, 
we evaluate the galaxy 3PCF at three indicative redshift 
values: 2 = 1,3 and 5 (for definiteness, we refer to them 
as "low" , "intermediate" and "high" redshift in what fol- 
lows). Specifically, this is assumed to be the redshift of 
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FIG. 2: Triangular configurations T4. The dashed (red) line 
shows a second triangle of the same class. 



the galaxy at location xi , while without loss of generality 
we always assume that the distances of the three points 
from the observer obey the condition Xi ^ X2 ^ Xs- 
Small corrections due to the different values of the growth 
factor at the three locations are easily included since they 
correspond to multiplicative factors at tree-level in FT. 
We compute the correlation function ^ and the auxiliary 
functions r\ and e from the power spectrum as described 
in the Appendix, assuming the matter transfer function 
provided by the CAMB code [51]. Therefore, acoustic 
features in the power spectrum are included. However, 
we do not include any smoothing effects due to non-linear 
gravitational clustering. These are beyond the scope of 
this paper, and should in any case be quite small at the 
redshifts considered. 

The last ingredient required to evaluate a realistic as- 
sessment of the lensing contribution to the measured 
galaxy 3PCF is a reasonable estimate of the bias param- 
eters 61, 62 and the number count slope s. We adopt the 
expressions for the bias parameters derived in the frame- 
work of the Halo Model and of the Halo Occupation Dis- 
tribution (HOD), [52, 53]. We will refer in particular to 
the values computed in [34] for several forthcoming high- 
redshift spectroscopic surveys, where the HOD assumed 
is the one proposed by [54]. They are derived in terms 
of the survey mean redshift and expected galaxy density; 
we assume rig ~ 5 x lO^"' h? Mpc^'^ for the z = \ and 
z = 3 examples, and rig ~ 50 x 10~^ Mpc~'^ for the 
z = 5 case. We will also consider values for the sur- 
vey volume oi V ^ Ah'^Gpc^, V = 2.7ft.~^Gpc^, and 
V — 3.4 Gpc'^ respectively for low, intermediate and 
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FIG. 3; Triangular configuration Tl. The upper panels show the absolute value of the unlensed 3PCF, the GGG group (black, 
long-dashed line), the contributions due to magnification, sum of the GGL, GLL and LLL terms (red, short-dashed line) and 
the total (blue, continuous line). The lower panels show the ratio between the magnification contributions and the unlensed 
3PCF (continuous line), again as absolute value. We also distinguish the GGL (short-dashed line), GLL (long-dashed line) and 
LLL terms (dotted line). First, second and third columns correspond respectively to low redshift {z = 1), intermediate redshift 
(z = 3), and high redshift (z = 5). Specific values for the bias parameters are included as discussed in the text. 



high redshift cases, which will be used later for an esti- 
mate of the expected error on the 3PCF measurement. 
These numbers roughly correspond to the characteristics 
of the proposed surveys WFMOS (0.5 < z < 1.3), HET- 
DEX (2 < z < 4), and CIP (3.5 < z < 6.5). Under these 
assumptions we obtain bias parameter values of bi 2± 2.0 
and 62 — 1-0 at z = 1, &i ~ 3.6 and &2 — 6.2 at z = 3 
and 61 2± 3.7 and &2 — 6.8 at z = 5. 

The value of the number count slope s strongly de- 
pends on the particular sample, i.e. type of object con- 
sidered as well as the magnitude cut applied. For galax- 
ies, typical values are in the range 0.2 < s < 0.6 for 
faint magnitude cuts, and 0.4 ^ s ^ 1.5 for bright mag- 
nitudes [44, 55]. For QSO, number count slopes range 
from s — —0.2 to s = 0.8 [56, 57]. In contrast to the 
bias parameters, the number count slope can be inferred 
from the observed distribution of galaxy apparent mag- 
nitudes. Furthermore, by changing the magnitude cut, 
it is possible to vary s within a limited range [57]. As 
a fiducial value, we assume s = 0.6 here, so that ci = 1 



and C2 = 2. 

In Fig. 3 (upper panels) we plot the absolute value 
of the different contributions to the 3PCF for z = 1, 
3 and 5 (left to right column) for the Tl class of tri- 
angular configurations. The long-dashed (black) curve 
represents the unlensed 3PCF, including the non-linear 
bias contribution, as a function of r, i.e. the line-of-sight 
separation between the points Xi and X3 (see Fig. 1, left 
panel), while the short-dashed (red) and the continuous 
(blue) curves represent the magnification contribution 
(all terms) and the total, lensed 3PCF, respectively. Note 
that, for this specific configuration, the 3PCF changes 
sign at about r ^ 30/i^^Mpc and several times in the 
range of scales corresponding to the acoustic features. 
This is due to the fact that the auxiliary functions rj and 
e in the tree-level PT expression for the 3PCF are related 
to derivatives of the matter correlation function which 
exhibits a local minimum and maximum in that range. 
In the lower panels we plot the ratio of the magnifica- 
tion contributions to the un-lensed 3PCF including the 
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FIG. 4: Triangular configuration T2, z = \, 3, and 5. 



bias. We also distinguish the contributions of the differ- 
ent groups identified in the previous section, that is the 
GGL terms (short-dashed curve), the GLL terms (long- 
dashed), and the LLL terms (dotted), while the solid line 
represents the sum of all of them. 

For most of the configurations and rcdshifts consid- 
ered, the magnification corrections to the 3PCF are above 
the percent level, especially at line-of-sight separations 
greater than ~ 50/i~^Mpc. Magnification corrections 
exceed 10% for separations above 40 Mpc at redshift 
z = 3, and above 30 Mpc at redshift z ~ h. As 
mentioned before, these trends are expected: while the 
magnification terms stay roughly constant, the unlensed 
3PCF is oscillating and decreasing in amplitude. This 
fact implies that the magnification corrections become 
more important for increasing separations, and that in 
the regions where the unlensed 3PCF changes sign the 
observed 3PCF will be completely dominated by lens- 
ing effects. Furthermore, it is apparent that the mag- 
nification effect increases with redshift. This fact has a 
straightforward physical explanation: the further away 
the sources are located, the more cosmological structure 
responsible for lensing is present between the source and 
the observer, which affects several of the correction terms. 
Finally, it is also clear from Fig. 3 that while the terms 



belonging to the GLL group, especially GLL-C, are dom- 
inant, the terms belonging to the LLL group have the 
strongest redshift evolution, and will eventually overtake 
the other terms at sufficiently high redshift. 

In Fig. 4 and 5 the same plots are presented for the 
triangular configurations T2 and T3, respectively. Most 
of the comments made about Tl are also valid for the tri- 
angular configurations T2 and T3. At intermediate and 
high redshifts, the magnification corrections exceed 10% 
for separations of 40 Mpc and eventually dominate 
the observed 3PCF. For T2 and T3, the unlensed 3PCF 
is again oscillating and decreasing in amplitude while the 
magnification corrections stay roughly constant. 

In Fig. 6 we finally present the results for the equi- 
lateral configuration T4. In this case, the magnification 
corrections behave somewhat differently as function of 
the separation. This is because in this particular case 
the transverse separation between the sources is not held 
fixed. Hence, the magnification correction is decreasing 
with increasing separation roughly in the same way as 
the overall 3PCF. In this case, then, the lensing contri- 
butions to the 3PCF become overwhelming only when 
the unlensed 3PCF is changing sign and therefore going 
through zero. As before, the lensing contributions in- 
crease with redshift: at r|| > 40/i~^Mpc, they amount 
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FIG. 5: Triangular configuration T3, z = 1, 3, and 5. 



to less than 1% for z = 1, a few percent at z = 3 and 
about 10% at z — 5. 

Going back to equations (18)-(31), it is clear that the 
magnification terms in different groups depend differently 
on the value of s (via ci and C2): the GGL terms scale as 
ci and C2, the GLL terms as and C1C2 while the LLL 
terms as cf and cfc2. Hence, depending on the value of 
s of the sample, different groups of terms will dominate. 
This is shown in Fig. 7, where we show the three sets of 
contributions as a function of s, for a fixed triangle at 
z = 3 (the Tl example with r = 60 /i"^ Mpc). Note that 
the GGL and LLL terms are odd with respect to (s — 0.4), 
i.e., they change sign for s < 0.4, while the GLL terms 
are even. The dominant individual contribution in this 
case turns out to be, as expected, the {5i52K2K3) (GLL- 
C) term, for all values of s. Note that this term increases 
rapidly with s around the value s ~ 0.6 chosen for the 
plots of Fig. 3-6. Thus, for a sample with high number 
count slope, the prospects of detection are significantly 
more optimistic than shown in Fig. 3-6. In the event 
of a detection, the scaling with s should be observed by 
measuring the 3PCF with different apparent magnitude 
cuts in the survey, corresponding to different values of s. 
In principle, this would allow to disentangle contributions 
with different scaling behavior, e.g. GLL and LLL terms. 



B. Reduced three-point function 

The reduced three-point function is commonly used to 
measure the non-Gaussianity of the galaxy distribution, 
as it contains all the information of the full 3PCF C, but 
is, to first order, independent of the growth and ampli- 
tude of matter fluctuations. It is given by 



(5(X1,X2,X3) = 



C(xi,X2,X3) 



C12C13 + ^12^23 + ^13^23 



(64) 



The reduced bispectrum is defined in an analogous way. 
In order to exemplify the effects of weak lensing on the 
reduced 3PCF, we show the intrinsic (unlensed) as well 
as lensed Q as a function oi 9 = <(ri2,ri3), the angle at 
vertex 1, for a triangle at z = 3 (Fig. 8). The lensed Q is 
shown for different orientations of the triangle, specified 
by the angle (3 between the normal to the plane of the 
triangle and the line of sight. /3 is chosen so that if /? > 0, 
vertex 3 is further away than 1 and 2, while it is closer for 
/? < 0. We choose ri2 to be perpendicular to the line of 
sight with ri2 = 20 h'^ Mpc, and ria = 80 h'^ Mpc. For 
this figure, we chose a value of s = 1.2 for the number 
count slope, corresponding to (5s — 2)/6 « 1.3. This 
value, while high, is not unexpected for bright magnitude 
cuts (e.g.. Fig. 1 in [44]). 
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FIG. 6: Triangular configuration T4, z — 1, 3 and 5. 



Following the discussion earlier, we expect the largest 
lensing contributions when /? = ±7r/2, i.e. when the tri- 
angle is oriented along the plane containing the line of 
sight, since in that case the projected separations are the 
smallest. In calculating the lensing corrections, we ex- 
cluded configurations where the perpendicular separation 
of any two points became less than 10 Mpc, since in 
those cases the perturbation theory approach is expected 
to break down. This results in the gap in the curve for 
/? = 7r/2. In that orientation ri-^± — > for 9 = tt/2. As 
the correlation function rises for r ^ 0, we expect the 
lensing contribution to rise sharply in this range. 

Apparently, the lensing effect on Q is appreciable in 
this configuration. Even for (3 = 0, i.e. a triangle lying 
in the sky plane, there is still a finite lensing contribution 
from the LLL terms. For inclined orientations, GGL and 
GLL terms contribute significantly. There is an enhance- 
ment of Q around 6 ^ tt/2 where vertex 3 passes behind 
(/3 > 0) or in front of {(3 < 0) vertices 1 and 2. For /3 
approaching +7r/2, there is noticeable structure in the 
lensing contribution. This is due to the sum over GGL 
and GLL terms peaking at different values of 9. 

Note that the observed Q will be given by the 3PCF 
including lensing effects divided by the product of the 
observed two-point correlation functions which will also 



include contributions from magnification bias. Since we 
focus on the lensing effects on the 3PCF, we did not 
include them here. 



C. Observability of the magnification contribution 

In order to assess the observability of the lensing effects 
discussed here, we have to estimate the expected statisti- 
cal uncertainty on the observed three-point function. In 
Appendix D, we derive a simple and very approximate 
expression for this uncertainty, assuming correlations are 
small (as is the case on large scales). To measure the 
three-point function for a fixed triangle, one defines bins 
in configuration space, so that for a given galaxy at ver- 
tex 1, triangles are counted for galaxies in the volume dV2 
around vertex 2 and dV^ around vertex 3. For a survey 
with Ng galaxies with a galaxy volume density of ng , we 
then obtain: 

a{0 = Ng^'\ I +36. (65) 
" y nj^dV2dV3 

For the rough estimates desired here, let us assume a 
typical bin size of dV2 dV^ 10* Mpc'^, corre- 
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FIG. 7: Ratio (absolute value, as in the lower panels of 
Figs. 3-6) of the lensing contributions GGL, GLL, and LLL to 
the unlensed GGG terms, as a function of the number count 
slope s. The triangle, Tl with r — 60 Mpc at z — 3, is 
kept fixed. 



FIG. 8: The reduced three-point function Q as a function of 
9 — <(ri2,ri3) for a triangle with ri2 = 20 Mpc, ns = 
80 h^^ Mpc at z = 3. The solid (black) curve shows the 
intrinsic Q, while the other curves show the total Q including 
magnification effects for different orientations of the triangle 
(see text). We assumed s — 1.2 for this figure. 



spending to a linear dimension of ~ 20 h^^ Mpc. Us- 
ing the values for the galaxy density and volume intro- 
duced above for the three ideal surveys at redshifts 2 = 1, 
z = 3 and z = 5 we obtain uncertainties on the specific 
configurations given, respectively, by a{C) ^ 4 • 10~^, 
cr(C) ~ 5 • 10"^ and cr(C) ~ 1.4 • 10"^ Note that these 
uncertainties correspond to triangular configuration bins 
with a given orientation, so that they can be directly 
related to the absolute values for the 3PCF plotted in 
Fig. 3, 4, 5 and 6. This suggests that magnification ef- 
fects should be potentially detectable in the galaxy 3PCF 
at redshifts below z w 3, and clearly measurable at z > 4. 

We expect, on the other hand, that when the correla- 
tion function is measured by averaging over all possible 
orientations for a given triangular shape, the magnifica- 
tion effect will be significantly reduced. However, in this 
case the expected error will be reduced as well due to 
the higher statistics. A proper assessment of the magni- 
tude of the correction given a specific 3PCF estimator is 
however beyond the scope of this work. 

Finally we should mention that another possibility to 
detect the lensing effects in much smaller surveys is given 
by cross-correlating galaxies or quasars at different red- 
shifts, or even across different surveys. In this case, there 
is no intrinsic correlation due to the large line-of-sight 
separation, while the magnification terms are still signif- 
icant at small transverse separations. Such a three-point 
magnification has indeed already been detected in the 
SDSS [56] in terms of the galaxy-quasar skewness. How- 



ever, no direct comparison between our results and those 
measurements can be made as we limit ourselves to the 
correlation function on much larger scales, where pertur- 
bation theory is applicable. 

IV. CONCLUSIONS 

In this work we analyzed how weak gravitational lens- 
ing affects the observed three-point correlation function 
of galaxies and quasars. As a congruence of photons trav- 
els from the source to the observer, the intervening dis- 
tribution of matter acts by focusing or defocusing it, thus 
altering the magnitude of the object observed. As a con- 
sequence of this effect, any magnitude limited sample - 
and any correlation function measured from it - will be 
somewhat contaminated by magnification bias. 

In a consistent second-order treatment in perturbation 
theory, weak gravitational lensing contributes 14 differ- 
ent terms. The lensing contribution to the 3PCF strongly 
depends on the configuration of the triangle considered, 
and in particular on its projection on the sky plane. If 
the projection is held fixed while the triangle is stretched, 
magnification bias can account for an appreciable fraction 
of the observed 3PCF at intermediate to high redshifts. 
On the other hand, for triangles oriented so that the pro- 
jection on the sky plane is of similar size as the trian- 
gle itself, the magnification corrections to the 3PCF are 
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smaller, reaching percent level at intermediate redshifts 
and of order 10% at high redshifts. The lensing con- 
tributions are a rapidly varying function of the number 
count slope s which depends on the sample considered. 
While we adopted an average value of s = 0.6 here, the 
lensing effects increase significantly for larger values of s 
(e.g.. Fig. 7 and 8). A rough estimate of the signal-to- 
noise suggests that the effect of magnification bias on the 
galaxy and quasar 3PCF may be detectable by itself at 
redshift 2 « 3 and above. 

The results presented here suggest that such magnifi- 
cation corrections need to be taken into account when 
measurements of non-Gaussianity are carried out using 
magnitude limited samples from galaxy and quasar cat- 
alogs at high redshifts. If neglected, they could in fact 
bias the determination of the bias parameters bi and 62, 
and affect the sensitivity of galaxy surveys to primordial 
non-Gaussianity. However, one can also argue that the 
distinctive dependence of the magnification effects on s 
as well as bi , 62 may be used in the future to unambigu- 
ously identify the lensing contribution, and to disentangle 
the different terms. Ultimately, it might be possible to 
use the lensing of the 3PCF as a tool to measure the 
intervening distribution of matter and the linear growth 
factor. A more detailed analysis of such aspects will be 
the focus of forthcoming work. 

Since this work represents just a first step in the in- 
vestigation and the possible exploitation of the lensing of 
the 3PCF. here we only considered a few specific triangu- 
lar configurations. However, one would expect that the 
lensing effects on the projected 3PCF, i.e. the three-point 
correlation of galaxies in a given redshift bin for a certain 
projected triangle, will be much larger. This is because 
many triangles with large line-of-sight separations within 
the redshift slice will contribute. On the other hand, by 
counting galaxies in a wide redshift bin, a spectroscopic 
redshift for each galaxy is no longer necessary, and much 
larger galaxy statistics become available. A generaliza- 
tion of this is to cross-correlate galaxies from different 
redshift slices, where the intrinsic correlation is negligi- 
ble. The lensing effects on the projected 3PCF will be 
studied in a future paper. 

Also, it would be desirable to extend the modeling of 
lensing effects to smaller separations, were they arc ex- 



pected and observed to grow rapidly [56]. In principle, 
one can use non-linear models for the matter power spec- 
trum in the expressions derived here [equations (37)-(43)] 
to extend the range of validity. However, this is not com- 
pletely consistent, since the separation into terms such as 
(Skk) , {SSkk) is in itself perturbative to leading order in 
the matter ovcrdcnsity. This caveat is equally valid for 
the calculation of magnification effects on the two-point 
correlation function [43, 44, 47]. A full non- linear calcu- 
lation of lensing effects on galaxy correlation functions 
would be considerably more involved and has not been 
attempted yet. 

Finally, the treatment presented here for the magnifi- 
cation bias of the galaxy and quasar 3PCF is a gener- 
alization of the one given in [43, 44. 58] for the 2PCF. 
We point out that the perturbation framework used in 
Sec. II B to derive the magnification corrections is com- 
pletely general: given enough patience and time, it is 
straightforward to extend the treatment given here to 
correlation functions of higher order. One might expect 
that higher order correlation functions will be more and 
more affected by magnification bias. 
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APPENDIX A: MATTER AND GALAXY 3PCF 
1. Conventions 

In all derivations, we use the following convention for the Fourier transform: 

<5k = J (^a^^'^'^W' '^W = / d^ke-'^^ ^'Sk. (Al) 
The power spectrum is defined as 

((5ki<5k.) =<5z3(ki+k2)P(fci), (A2) 

while the bispectrum is defined as 

{Si.,6k,6^,) = SoCki + k2 + k3)B(fci, ^2, h). (A3) 

The dependence of P(A:i) and of B{ki, fc2, fcs) on the redshifts/comoving distances is implicitly understood and omitted 
throughout the text. Since we are consistently using the linear power spectrum, such dependence can simply be 
accounted for by multiplying by the appropriate growth factors. 

To obtain expressions in the other commonly used Fourier convention, where (5(x) = J d"^fc/(27r)'^(5k exp(zk • x), just 
replace P{k) with P{k) / . The same is true for the correlation functions defined below. Further, we again use 
the notation Xy = x^ — x^-, Xij = |xi — Xj|. 



2. Computation of the matter 3PCF — Auxiliary Functions 

The matter two-point correlation function is obtained from the matter power spectrum as 

i{x) =4n [dk fc2p(fc)!i!ilM. (A4) 
J kx 

Assuming the tree- level expression for the bispectrum given in equation (1), Jing and Borner [48] showed that one 
can write 

C(xi, X2, X3) = yC(a;i2)C(a;i3) - [?72(a;i2)??o(a;i3) + '7o(a;i2)'72(a;i3)] X12 ■ X13 

+ y [e(a;i2)e(a;i3)(xi2 ■ Xi3)^ + ry2(a;i3)e(a;i2)a;i2 + 'n2{xi2)e-{xi:i)x\^ + 3r/2(a;i2)r/2(a;i3)] 

+2 perm., (A5) 

where 

/ N . f , P(k) kx cos(kx) — sm(kx) , . 

rj^ix) ^ A^j dke^ "^-^^ L^, (A6) 

e{x) ^ in [ rffcP(fc)3[^i^(fc^) - kxcosjkx)] - fc^x^ sin(fcx) _ (^7) 
J k X 

This expression is quite convenient since it allows us to determine the 3PCF from the functions ^, rji and e which have 
to be computed only once. These functions are plotted in Fig. 9. 



3. Transverse 2PCF, 3PCF and auxiliary functions 

When applying the Limber approximation, one can obtain very similar expressions for the magnification con- 
tributions in terms of the projected two- and three-point correlation functions, ^ and Cj respecively. The projected 
two-point function ^ is defined as the 2-dimensional Fourier transform of the power spectrum in terms of the transverse 
wavenumber kj^ (equivalent to forcing kz =0): 

C{x^) EE (27r) J d^k^P{k^)e-'^^-''^ (A8) 
= {2TTf I dkkP{k)Jo{kx±), (A9) 
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FIG. 9: Linear correlation function ^(a;) and auxiliary functions rjo{x), rj2{x) and e{x) as a function of x, evaluated at redshift 
z — 0. In the left panel, the non-linear correlation function is also shown. The right panel shows the absolute value of the 
functions. 



where Joix) is the 0th order Bessel function. The transverse three-point function ^ is defined analogously as 

C(xii,X2±,X3±) = {2TTf Jd^h^e-^'^'^-^'^ Jd^k2±e-'''-^-^'^B{ki^,k2±)- (AlO) 

Note that ^ has dimensions of length, and C of length squared. In order-of-magnitude terms, ^{x±) ~ x±£^{x±), 
Cix±,x±) ~ x']_({x±, x±). In analogy with what is done above for the 3-dimensional 3PCF, one can write 

C(xii,X2±,X3i) = -^C{xi2±)^ixi3±) - [rj2{xi2±)rjo{xi3±) +rjo{xi2±)rj2{xi3±)] Xi2_L • Xi3_l 



+ y [e(a;i2j_)e(a;i3J_)(xi2± • xi3x)^ + ri2{xi3±)e{xi2±)x 
+V2ixi2±)e{xi3±)xlrij_ + 2rj2{xi2±)T]2ixi3±)] + 2 perm. 



2 

12± 



where 



mix) 



fc' 



Ji{kx), 



(All) 

(A12) 
(A13) 



Here, Ji{x) denote Bessel functions of integer order. These functions arc plotted in Fig. 10. Since all auxiliary 
functions are proportional to the linear matter power spectrum, they scale with the square of the matter growth 
factor. For brevity, we suppress the redshift- or ^-dependence of these functions here and in the following. 



I{x) = I dkP{k)[kxJo{kx) ~2Ji{kx)]. 



APPENDIX B: DERIVATION OF CONVERGENCE AND SHEAR CORRELATORS 



This appendix contains the derivation of convergence k and shear 7 as well as representative derivations of some of 
the results shown in equations (37-43). 
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1. Convergence and shear 



The projected Icnsing potential ^ is defined as 

*(x) = - r dxWL{xs,x) mvix); x) - *(y(x); x)] , (Bi) 

Jo 

where $, ^E* are the gravitational potentials and \1/ = — $ at late times in ACDM and most dark energy models. Xs is 
the distance to the source, and y(x) = (x/xs)x parametrizes the line of sight. The convergence k(x) and the complex 
shear 7(x) are then given by 

2k = «',ii + «',22, (B2) 
27 = ^-.n - ^',22 + Zi^'.iz. (B3) 

We now use the Poisson equation. 

kH ^ ^ = - 5, (B4) 

z a a 

to write k and 7 in terms of the matter overdensity S. For the convergence, this yields: 

k{^) = C prfx-^^^^^^<5(y(x);x) (B5) 
Wl{Xs,x) I d'k^~^^yix)s^^^), (B6) 

Jo 

For the shear, we obtain: 

7(x) = C r dx^^%^ / d3fc^^^M±^,,(^)e--y(x). (B7) 
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2. Derivation of correlators 



As shown in section II, the different correlators involving k and 7 simplify considerably when using the Limber 
approximation. We now show this in more detail in a few sample derivations. First, consider ((5(xi)k(x2)), where we 
assume as usual xi < X2- 

(5(xi)^(x2)) = c r dx^^'^^^a\^i-y2\) 
c 

The first expression can be used to evaluate this correlator directly. However, in the fc-space formulation, we can apply 
the Limber approximation: the typical longitudinal separations X12 — y2z are much larger than the perpendicular 
separation, hence the exponential will oscillate rapidly unless < l/(xi2 — y2z)- Thus, we set k^ = 0, which yields 
a Dirac J- function in % — xi, fixing the lens to be at the same distance as point 1: 



a(x) 

Jo J 



a(Xi) J 



(B9) 



= c ^^^f'f ^^ g(xi2j.), (BIO) 
a(xi) 

where we have used the definition of the projected two-point function, equation (A8). 

In a similar way, we obtain expressions for the correlators {6"/) (again assuming xi < X2): 

(^(x,)7(x.)) = C r dx^^^^ I d^fcP(fc)^i^M±^e--(---). (BID 
Jo '^IXj J k 

Using the Limber approximation and setting fc^ — > k\ in the denominator as well, we find 

(<5(x07(x.)) 2.C^^^^^ [ d^fc^P(fc^)^L^4±^e->^-(-— ). (B12) 
a(Xi) J fcl 

The integral can be expressed as 



kl~kl+ 2ikik2 

k' 



2n I d^k^P{k S' '^^J-^'^^-^' e 



1. 

ikj_x± cos <! 



27r j dkx_k^P{kx_) j d0[2 cos^ - l]e" 

dk^P{k±_)[k^x^Jo[k^xi_) -2Ji{k^xi_)] = x\e{xi_), (B13) 



2ti 

where the function 't{x±_) is defined in equation (A13). Hence 



L A. ^Wl(X2,Xi) 2 



(J(xi)7(x2)) ^ C^^f^o: 2,i?(.Ti2,x). (B14) 



Similarly, for the (K7) term we get: 



(.(x07(x2)) =C^r dx^^^ r / d^kPik)'l^^±^e--iy^'y^\ (B15) 

Jo a(x) Jo a(x') J k^ 

which, applying the Limber approximation, yields 

/ / N / N\ L.A 2 , W^l(X1,X)W^l(X2,x) 2 ~f N /-niR^ 

(k(xi)7(x2)) ~ c / dx :rr-. yi2.±^iyi2,±)- (bi6) 



a2(x) 

In all cases, the auxiliary functions are to be evaluated at xi, or at x if inside an integral over x- Fig- H 
shows the four two-point correlators as function of the transverse separation xi2±, for Xi at z = 3, and keeping 
Axi2 = X2 ^ Xi fixed. Both the exact expressions and the Limber approximation results are shown (see the next 
section for a discussion) . In the case of (Sk) , (^7) , we also show the functions for different values of the longitudinal 
separation Axi2 = X2 — Xi- Clearly, these correlators increase strongly for increasing longitudinal separation, as they 
are proportional to Wl{x2tXi) ~ ^Xi2- For (kk) and (^7), the dependence on the longitudinal separation of 1 and 
2 is very small (as long as A.X12 ^ Xii X2), since it only enters the arguments of the lensing window functions in the 
integral. 
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FIG. 11: (5(xi)k(x2)), ((5(xi)7(x2)), {K(xi)fi:(x2)) and (fi:(xi)7(x2)) as function of the perpendicular separation xi2±. Here, 
xi is kept fixed at z = 3, while X2 is separated by Axi2 = X2 - XI = 200, 50, -50 h-^ Mpc. In the last case, X2 is closer 
than Xi. Shown are the results for the exact integral (in case of (^7) and (^7), the first-order calculation) and when using the 
Limber approximation (L.A.). 



APPENDIX C: LIMBER APPROXIMATION 



1. Discussion 



The Limber approximation is very useful to gain physical insight into the magnification effects, and for order-of- 
magnitude estimates. However, before using the simplified expressions for exact calculations, one should assess the 
accuracy of this approximation. We calculated the exact integrals for all of the terms in equations (37)-(43), except 
the (kkk), (<57), and (^7) terms. For the latter two terms involving the shear, the exact calculation is very difficult, 
and we include the first order correction which is good to 0{k^/k\_) instead. The calculation of this correction is 
outlined in the next section. Note that the first order correction to the shear terms are quite small, so that there is 
no need to go to higher orders. 

Fig. 11 shows a comparison of the exact expressions and the Limber approximation. For (i5k) and ((S7), the Limber 
approximation is not accurate anymore at large transverse separations. Since we keep the longitudinal separation 
Axi2 fixed, large Xi2± correspond to more and more transverse X12. Correspondingly, the deviations become smaller 
for larger Axi2, as is apparent from Fig. 11. The reason for this is that the Limber approximation assumes that 
separation vectors with < x± contribute to the correlation. However, if a;i2_L > Si2||, then this overcounts the 
contributing modes and overestimates {5k), {5j). As we lower Xi2\\ = Axi2 to and negative values, the result of the 
Limber approximation goes linearly to 0, while there is actually still a finite correlation at negative Axi2, i.e. when 
we correlate an overdensity at Xi behind the lens magnifying X2 . Assuming the lens is a distance A.xl in front of X2 , 
an order of magnitude estimate gives {S1K2) {HoAxL){HoXi2i_)i{rL), where = \/a;^2± + i^XL + |Axi2|)^- For 
xi2± = 100 h-^Mpc, Axi2 = -50 /i"^ Mpc, and setting Axl ~ 50 h'^Mpc, this yields \{6iK2) \ ~ 10"^. This is in 
rough agreement with the exact calculation, shown by the dash-dotted curve of {S1K2) for Axi2 = —50 in Fig. 11. 
Clearly, {S1K2) becomes very small for negative Axi2 and in fact decreases rapidly for increasingly negative values. 
Terms like this would only contribute for almost transverse separations, where the (Sk) term is in any case very small. 
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This justifies the neglection of terms involving (^2^1) etc. in the derivation of equations (18)-(31). 

For the (kk) and (k'j) terms, the Limber approximation is generally very good (we do not show them for other 
values of Axi2 as they are essentially the same). For (kk), there are very small deviations at small Xi2±, corresponding 
to almost longitudinal x. In this case, the Limber approximation underestimates the correlation, counting only purely 
transverse modes {yi2\\ < yi2±), while in the exact calculation non-transverse modes contribute as well. 

To summarize, we found that the Limber approximation holds well for quantities that, after the approximation, 
still involve an integral over the line of sight, like {kk), (^7), and {kkk). The integral serves to average out the effect 
of the approximation. For terms that only involve a single evaluation of the correlation functions, e.g., (Sk), (Sj), 
(Skk), the Limber approximation shows deviations, especially when the Limber-approximated expression goes to zero, 
while the actual correlation is finite. For these quantities, the exact integral should be used. 



2. Shear terms beyond the Limber approximation 



While the evaluation of the exact integrals is straightforward for the terms involving the convergence k, going 
beyond the Limber approximation requires more work for the shear terms - however, it provides some insight as well. 
Going back to equation (Bll) and using cylindrical coordinates {k±,kz,(j)), we can write {61^2) as: 

(<5(xi)7(x2)) - C £^ dx ^^^^g /(r ^ y2 - xi , x) (CI) 

/(r, x) = J k^dk^ J dkz P{k ^ ^kl + kl; x) J d^^l_^2+^^hj^ exp (-*(k^ -r^+k. r,, )) (C2) 

= J k±dk± J dfcz exp(— i/c^ r||) F(fc; x)^ J di^ (2 cos^ — 1 -f 2i cos 0sin 0) exp(— z(k^ • r^)) (C3) 

k^dk^ J dk, cxp(-zfc,r||)P(fc;x) ^2^|^2 (^^) [joik±r±) - -^Jiik±r^)^ , (C4) 

where we have aligned k_L so that (j) is the angle between k_L and r_L. The last integral is difficult to evaluate. In order 
to make progress, we use the fact that typically ry ^ r±, so that kz should be much less than k± to get a significant 
contribution. Thus, we expand the integrand in powers of fc^/fc^. For the power spectrum, this expansion yields: 

m = P(fc.) + ^^k^ofc.^ + o(^) . (C5) 

Here and in the following, we have suppressed the x-evolution of the power spectrum. Note that in this expansion 
there is no linear (and in general no odd) term in kz- We now obtain 



/(r) ~ / k^dk^ Jdkz exp(-ifc,r||) (^P{k^) + J-^l.^^^kl) (l - 
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x{2n) (^Jo{k±r^) - ^^J^{k^r^)^ . (C6) 

We then obtain zcroth and first order expressions for the integral, the first order expression being good to (fc^/fc^)*: 

/'"Hr) = j k^dk^P{k^)(2TT)(^J^{k^r^)- j^,h{k^rS^ J dkz cxp{~ikz ^) (C7) 

/«(r) ^ Jk.dk, + U^-^o)(2.)(jo(fc...)-^J.(fc...)) 

X J dkz exp(— ifc^ r||) (C8) 

The integrals over kz result in Dirac 5]j functions and derivatives. We sec that the zcroth order term is exactly what 
we obtained earlier in the Limber approximation, equation (B14): 

/(0)(i.) = IM! f dkP{k) [kr^ Mkr^) - 2Ji(ftr^)) = rl^r^) doi^), (C9) 
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where r\\ = X ~ Xi in equation (CI), so that the 5]j function fixes X = Xi- For the first correction to the Limber 
approximation, we get 

l''Hr) = / m I- -l) (kr^ Mkr^) 2J,{kr^)) ■ ^<5^(r||) ee V{r^) ■ -^5n{n\ (CIO) 

^ J \ / II II 

defining another auxihary function ^{r). By partial integration in equation (CI) for {5^) and equation (B15) for (^7), 
we can calculate the first correction to the Limber approximation for the shear terms, {5^)^^\ {kj)^^'>: 

/ \(i) ^2 , (Wl(x2,x)Wl{xi,x)\ ~, X .^1o^ 

where j/i2_l = {xlxi)^i2i^ and we have made the dependence on the evolution x explicit again. The derivatives of the 
lensing weight factors can be done easily, yielding polynomials in xi, x a-nd factors of H{z) and dH/dz. 

The efi'ect of the first order correction on (61^2) and (K172) is shown in Fig. 11 (right two plots). While for {61^2) 
there is some effect at large a;i2_L (similarly to {S1K2)), the correction to (K172) is negligible. In any case, we conclude 
that the first-order calculation is sufficient, and there is no need to go to higher orders in equation (C4). 



APPENDIX D: A ROUGH ERROR ESTIMATE FOR THE THREE-POINT CORRELATION FUNCTION 

In order to estimate the error on the three-point correlation function (3PCF), we start from the number of triples 
DDD observed, i.e. the number of triangles in a certain geometry bin formed by galaxies in a survey. Say the survey 
observes a total volume V containing Ng galaxies. We specify a triangle bin in such a way that, after fixing point 1, 
we count all triangles where a galaxy is found in the volume element dV2 around point 2 and a galaxy found in dV^ 
around point 3 (points 2 and 3 being determined by the triangle geometry considered). Then, the expectation value 
of triples is given by: 

(DDD) = 7V3 [1 + ^('^12) + ^('^13) + ^(^^23) + C('^i2, ri3, r23)], (Dl) 
N, ^ Ng {ngdV2) ingdVs) = N^^^, (D2) 

where Ug is the average volume density of galaxies, estimated by Ug — Ng/V . An estimator for the 3PCF C, is given 
in [59], § 55: 

; DDD -DDR „ 
RRR 

where D stands for points drawn from the data, while R stands for points drawn from a random distribution (with 
the same boundaries). Calculalting the variance of this estimator, we obtain: 

a^iO EE (C)2 (D4) 

{DDDf\ /DDD\^ ^ I DDDDDR\ 



[RRRf I \ RRR I \ {RRRY / 

I DDD \ I DDR \ l {DDRf \ I DDR \^ 

^ \RRRl\RRRl^\{RRRfl \ RRR / ' ^ ' 

To proceed further, we make the following approximations. The RRR in the denominator is assumed to be uncorrelated 
with the numerators. Further, we set: 

{{RRRf) = {RRR}'^ = Ni. (D6) 

In addition, if we assume that the galaxy correlations are weak, i.e. that (DDD) « N3, then the second and third set 
of terms in equation (D5) vanish. We then obtain: 

'^'(C) = 4 {{(DDDf) {DDDf) 
^^3 



aHiDDD)) 



Ni 



(D7) 
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In this approximation, the error on the 3PCF is simply the relative error on the number of triples DDD. The 
latter was calculated in [60] and is given to good accuracy by: 

<y^{{DDD)) = (DDD) + ^{DDD)\ (D8) 

The first term is the standard Poisson term expected for any number count. The second one was found for the first 
time in the derivation of [60]. Now we again neglect the galaxy correlations in equation (Dl) and set {DDD) k, Nt, = 
Ng ■n?gdV2dVz, to obtain: 



^(0 = JtTT^ + - N-'''\ hTl^^ 36. (D9) 



{DDD) Ng <^ y nldV2dVz 
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The first term again is the Poisson term: cr(C)poisson = ^triples' analogous to cr(^) « -^palrC '"^ found in [59], § 48, for 
the two-point correlation function. The second term in equation (D8) was added in [60], along with a similar term 
for the variance of pair counts. 



